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Lya Emitters at redshift z ~ 4.5 in the Extended CDF-S region: 

I, Lya Luminosity Function 
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ABSTRACT 

We present a spectroscopically confirmed Lyman alpha emitter sample at z ~ 
4.5 in the Extended Chandra Deep Field South region. Starting with spectra of 
68 candidate Lya emitters, we identified 46 z ~ 4.5 Lya emitters based on single- 
line detections with no continuum emission blueward of the line, and 5 Lyman 
break galaxies. The Lya confirmation rate is ~ 70% for candidates selected in 
our two narrowband filters, and the remaining candidates including 5 Lyman 
Break Galaxies and one strange multi emission line object. We did not find any 
significant CIV emission line or Hell emission line in the optical spectra even 
with the coadded optical spectra. A majority of the LAEs show asymmetric Lya 
line profile, and the coadded Lya line has wavelength-based asymmetry of a a 
= 1.491q^3 ^'^'^ flux-based asymmetry of a/ = 1.451q'q7. Deep X-ray, radio and 
the optical spectral data are also used to constrain the AGN contribution in our 
sample, only one detection (not targeted in spectroscopic observation) was found 
at z=4.5 in both X-ray and radio, while the other objects show non-detection even 
with a stacking method. Finally, the Lya luminosity function for our z ~ 4.5 
sample is consistent with others published in the literature for the redshift range 
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of 3 < 2; < 6, while the Lya luminosity functions in our two separate narrowband 
filters show a ~ 2(T difference. This difference is comparable to the differences 
among z = 5.7 and z = 6.5 Lya luminosity functions reported by Ouchi et al. 
(2010), which are thought to be an implication of cosmic re- ionization. 

Subject headings: galaxies: active — galaxies: high-redshift — galaxies:starburst 

1. INTRODUCTION 

Thanks to efficient wide field cameras, the number of star-forming galaxies known in 
the early universe has grown rapidly in the past decade. These high-redshift star-forming 
galaxies are selected mainly through two techniques, the dropout technique and Lya-line 
search technique. The former is known as the Lyman-break technique applied on the deep 
broadband images (e.g., Steidel et al. 1996), whereby high-redshift galaxies are identified 
via a flux discontinuity caused by Lyman limit absorption, and the latter is designed to 
search for the strong Lya line from the narrowband images as the Lya line redshifted to the 
windows of low night-sky emission (e.g., from z = 2.1 to z = 7.7, Cowie & Hu 1998, Guaita 
et al 2010, Gawiser et al. 2007, Rhoads et al. 2000, 2003, Ouchi et al. 2008, Wang et al. 
2005, Dawson et al. 2004,2007, lye et al. 2006, Finkelstein et al. 2008, 2009b,c and more). 
The Narrowband surveys of high-redshift Lya emitter(LAE) galaxies can reach a successful 
fraction of > 70% (Dawson et al. 2007, Wang et al. 2009). Lya searches at corresponding 
redshift can identify galaxies with strong Lyo; line emission, which might be too faint to be 
detected by the Lyman-break technique. 

Although Narrowband Surveys have found thousands of Lya emitters from z = 2.1-6.96, 
fewer than one thousand of them have been spectroscopically confirmed. It is very interesting 
that the Lya luminosity functions of LAEs do not evolve from redshift z~ 3 to z~6.5 (e.g., 
Malhotra & Rhoads 2004, Dawson et al. 2007, Wang et al. 2009, Ouchi et al. 2008, Ono 
et al. 20 11). Optical spectra of LAEs typically show asymmetric Lya line profiles (Rhoads 
et al. 2003). No CIV emission line (indicating of active galactic nuclei [AGN] activity) or 
Hell emission line (indicating of the first generation stars, as well as some kind of AGNs) 
was reported in the previous LAEs' spectra (Dawson et al. 2007, Wang et al. 2009). Deep 
X-ray images are the most effective ways to find AGNs. Previous 170ks and 180ks Chandra 
X-ray exposure in the Large Area Lyman Alpha (LALA, Rhoads et al. 2000) survey did 
not find any X-ray individual or average detections at z~4.5 (Malhotra et al. 2003, Wang 
et al. 2004). However, in the 250ks Extended Chandra Deep Field South (ECDFS) region, 
where exist the deepest X-ray field of 4Ms Chandra Deep Field South (CDF-S), we found 
one quasar from 112 LAE candidates at z~4.5 (Zheng et al. 2010). Excluding that LAE- 
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AGN, the remaining did not show any detection even with stacking methods. At redshift 4 
< z < 2, more AGNs are found through deep X-ray surveys, but the AGN fraction is still 
<5% (Gawiser et al. 2007, Guaita et al. 2010, Ouchi et al. 2008). This is quite different 
compared with the z ~ 0.3 LAEs detected through GALEX, where the AGN fraction is as 
high as 15%-40% (Finkelstein et al. 2009a, Scarlata et al. 2009, Cowie et al. 2010) 

Here we report spectroscopic observations of z ~ 4.5 Lya emitters in the 0.34 deg^ 
ECDFS region. This region has an extraordinary amount of complementary data, including 
high-resolution optical images as the Great Observatories Origins Deep Survey (GOODS; 
Giavalisco et al. 2004) from the Hubble Space Telescope, deep ground-based UBVRIzJHK 
photometry from MUSYC and EIS, mid- and far- IR observations from Spitzer, GOODS, 
EIS and MUSYC, deep radio data from VIA array (Miller et al. 2008, 2011), and deep 
X-ray data from Chandra (Giacconi et al. 2002; Alexander et al. 2003; Lehmer et al. 2005; 
Luo et al. 2008; Xue et al. 2011). Our candidates and spectroscopic targets are mainly from 
two narrowband images designed for z ~ 4.5 LAE hunting by Finkelstein et al. (2009c). We 
present our photometric and spectroscopic observations in §2, spectroscopic results in §3, 
and discuss the AGN contamination fraction, the constraints on Population III stars, and 
the properties of luminosity functions in §4. This paper (paper I) contains the main results 
of the spectroscopic surveys. Throughout this work, we assume a cosmology with Hq = 71 
km s"^ Mpc"^ (Wright et al. 2006), = 0.27 and = 0.73. At redshift z = 4.5, the age 
of the universe was 1.36 Gyr, with a scale of 6.7 kpc/" and redshift bin of 6z = 0.03 implies 
a comoving distance of 18.75 Mpc. Magnitudes are given in the AB system. 



2. OBSERVATION 

2.1. Candidate Selection through Narrowband Imaging 

The LAE candidates were selected using narrowband imaging of the GOODS Chandra 
Deep Field South (CDF-S; RA 03:31:54.02, Dec -27:48:31.5, J2000) obtained at the Blanco 
4m telescope at Cerro Tololo InterAmerican Observatory (CTIO) with the MOSAIC II cam- 
era. Three 80 A wide narrowband filters were utilized to obtain deep narrowband images 
(2.75 hr of observation in NB656, Finkelstein et al. 2008; 5 hr in NB665 and 5.5 hr in NB673, 
Finkelstein et al. 2009c). The LAE candidates are selected in Finkelstein et al. (2009c), 
based on a 5 a significance detection in the narrowband, a 4 a significance narrowband flux 
excess over the R band, a factor of 2 ratio of narrowband flux to broadband flux density, and 
no more than 2 a significant flux in the B-band. The first three criteria ensure a significant 
line detection, while the last criterion checks that it is at z > 4. With these creteria, we have 
selected 4 candidates in NB656 (Finkelstein et al. 2008), 42 in NB665 and 85 in NB673 (11 
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candidates were selected in both NB665 and NB673, Finkelstein et al. 2009c)lj. Candidates 
with GOODS B-band coverage were further examined in GOODS B-band image, and 3 can- 
didates in NB665 and 5 candidates in NB673 with significant GOODS B-band detections 
were excluded (see the tabled]). So finally we select 112 LAE candidates at z ~ 4.5. 



2.2. Spectroscopic Observations 

Our spectroscopic data were obtained using the IMACS (Dressier et al. 2006) short 
camera (//2, with a 27'. 2 diameter field of view) on the 6.5 m Magellan I Baade Telescope in 
2009 September 10-ll(through Steward Observatory time, PI Rhoads) with the 200 line/mm 
grism and 2009 November 11-12 (NOAO PID 2009B-0371, PI Finkelstein) with the 300 
line/mm grism. The 200 line/mm grism has a Xuaze = 6600 Aand a resolution of 2.037 A 
pixel"^ with a range of 4000-10500 A, and the 300 line/mm grism has a same Xuaze = 6700 
A, and a resolution of 1.341 A pixel"^ with a range of 4000-9200 A. Five multislit masks 
(see Table 2) were observed for 2.8-4.5 hr in 0.5 hr increments. The masks have slit widths 
of 0.8 arcsec. Each slit mask included approximately 20 candidate Lya emitters (mixed 
in with roughly ~150 (with 200 line/mm grism) or ~50 (with 300 line/mm grism) other 
spectroscopic targets). Of these, 16 candidates were covered by more than one mask. In 
total we had targeted 64 LAE candidates. 

The data were reduced using the IMACS version of the Carnegie Observatories System 
for MultiObject Spectroscopy (COSMOS) data reduction package!. We first determined 
two-dimentional wavelength solutions for each science exposure using arc lamp exposures 
taken immediately before or after. The wavelength residuals in the calibration is about ~2 
pixels. 

After wavelength calibration, each frame was first bias-subtracted and fiat-fielded. We 
then performed sky subtraction following the algorithm described by Kelson (2003), and 
extracted one-dimensional spectra from the two-dimensional spectra using the task "extract- 
2dspec" for each slit. 

To control for possible spatial shifts along the slits between individual exposures of a 



^Thc NB656 candidates were selected only in the overlap area between the MOSAIC image and the 
GOODS CDF-S data, which is why only four objects were selected. The other two catalogs consist of all 
selected candidates over the overlap region between the MOSAIC image and the ESO Imaging Survey, which 
consists of a much larger area. 

^http://obs. carnegiescience.edu/Code/cosmos/Cookbook. html 
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Table 1. Number of Candidate, Targeted and Confirmed Lya Emitters at z ~ 4.5 in Each 

Narrow Band Filter in the CDF-S field. 



Filter 


Seeing 




Candidates'' 


Targeted'^ 


Confirmed'^ 


NB656 (HaO) 


0".90 


24.9 


4 


3 


3 


NB665 (Ha8) 


0".92 


25.0 


33 


17 


11 


NB673 (Hal6) 


0".91 


25.2 


75 


44 


32 



^5-0" limiting magnitude at corresponding narrowband. 

'^11 LAE candidates were covered by both Ha8 and Hal6 filter, and we 
had divided them into corresponding narrowband (5 into NB665 candidates 
and 6 into NB673 candidates) due to their narrowband flux. 

'^G of 11 LAE candidates covered by both Ha8 and Hal6 filter were targeted, 
3 were confirmed as NB665 LAEs and 3 were confirmed as NB673 LAEs, and 
the confirmation is consistent with the narrowband division. 



Table 2. Number of Targeted and Confirmed Lya Emitters at z ~ 4.5 in Each Mask in 

the CDF-S field. 



Mask Grism (line/mm) Exp. Time (ks) Targeted Confirmed 



1 200 10.2 20 11 

2 200 12.2 26 20 

3 300 16.2 22 15 

4 300 14.4 16 15 

5 300 14.4 10 7 
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mask, we measured the trace locations of the brightest continuum sources in the mask. We 
corrected for any measured shifts while stacking the exposures for each mask to increase the 
quality of the stacked two-dimensional spectra. We also identified and removed cosmic ray 
hits while stacking the multiple exposures of each mask. 

3. SPECTROSCOPIC RESULTS 

3.1. Line Identification 

We identified 46 single emission lines as Lya-emitting galaxies at z ~ 4.5, among our 
64 LAE candidates. In the remaining 18 candidates, 1 is possible (NB673-27) due to its 
strange spectra, 3 lost by CCD gap (NB665-15) or contaminated by skylines (NB665-32 and 
NB665-38) from nearby slits, 9 are undetected, and 5 show continuum breaks with little- 
to-no Lya line. In this work, the fraction of spectroscopic confirmation is about 70%-80%, 
and the contamination fraction is about 14%-19%. We present a catalog of the 46 confirmed 
Lya-emitting galaxies at z ~ 4.5 in CDF-S region in Table |3l and their redshift distribution 
is plotted in figure El There were 11 candidates detected in both NB665 and NB673 images, 
and 6 of 11 were targeted, and we had divided them into corresponding narrowband due to 
their narrowband flux. 3 were confirmed as NB665 LAEs and 3 were confirmed as NB673 
LAEs, and the confirmation is consistent with the division (see figure |3]). In Figured! we 
present narrowband stamps of all the candidates, targets and confirmed LAEs. Their Lya 
flux distributions are plotted in Figure |2J It is well seen that all the targets with f^j^ct > 
3.7x10"^'^ ergscm~^s~^ are confirmed as LAEs. 

Beyond the multiple emission lines in the low-redshift interlopers, we may identify Lya 
emission by its characteristically asymmetric morphology. Following Rhoads et al. (2003, 
2004), we perform two measurements of the line asymmetry: the wavelength-based asym- 
metry defined as ax = {\io,r — \)/{.\ ~ ^io,b)^ and the flux-based asymmetry defined as 
a/ = f\d\)/{j^^^^f\d\), where the Ap is the wavelength of the emission-line peak, 

and the Aio,b and Aio,r are the wavelengths where the flux density equals 10% of the peak 
on the blue side and red side, respectively. In figure IH we plot a\ versus aj for the 46 
confirmed LAEs. The error bars on a\ and are estimated based on a large number of 
Monte-Carlo simulations, in which we added random noise to each data bin. There are 34 
LAEs show both a\ and aj greater than 1, and only 2 LAEs show both a\< 1 and a/< 1 
even with their 1-cr error. These two LAEs (NB665-14 and NB673-72) have extended 2-d 
spectra, so the asymmetry of their Lya lines are hard to measure. We also notice that the 
asymmetry measured with the 300 line/mm grism show less dispersion than that measured 
with 200 line/mm grism, this is mainly due to the more exposure time observed with 300 
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Fig. 1. — Narrowband stamps (size: 10"x 10", 0.27"per pixel) of all LAE candidates, 
including confirmed LAEs with red circles, un-confirmed candidates with blue circles, and 
un-targeted candidates without any marks. There are 7 objects marked as blue circles 
without red-dashed diagonal lines, except 673-27 which shows strange profile and is hard to 
judge, the other 6 objects show continuum with little-to-no Lya emission line. All stamps 
are divided into 13 lines (line 2 is empty, and line 11, 12 have 27, 26 candidates, respectively), 
with a decrease bin of Algio(/Lya) of 0.07 dex for nearby two lines. It is well seen that the 
candidates with larger fLya are more likely LAEs at z ~ 4.5. 
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fL,„ (10-1' erg cm-^) 



Fig. 2. — Lya flux distributions of all candidates (empty histogram), IMACS targets (grey 
filled histogram) and confirmed Lya emitters (blacked filled histogram) . Candidates selected 
in different narrowbands (NB665 and NB673) are plotted in the middle and top panels. 
The 11 LAE candidates covered in both NB665 and NB673 are divided into corresponding 
narrowband based on their Lya fluxes. 
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Fig. 3. — Narrowband filter transmission curves (low-right axis) and redshift distributions 
of confirmed LAEs (top-left axis). Redshift distributions confirmed with filters NB656, 
NB665, and NB673 (from left to right) show same colors as the filters' transmission curves, 
Orange and light green filed histograms are confirmed candidates observed in both NB673 
and NB665 filters. The NB656 candidates are selected with a much smaller area and shallower 
narrowband image compared with NB665 and NB673, 
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line/mm grism than with 200 hne/mm grism, also due to the better spectral resolution for 
300 line/mm grism than that for 200 line/mm grism. 

3.2. Spectroscopic Calibration 

We use one standard star UID1147 (Pirzkal et al. 2005) to calibrate the spectra. The 
comparison of Lja line fluxes integrated from spectroscopy and estimated from photometry 
is plotted in the figure |5l A large dispersion is seen in the fiux-fiux plot, and line flux from 
optical spectra with 10.2 ks exposure of grism-200 are systematic larger than that from 
photometric data, which might be due to the slit-loss of the standard star observation. In 
the following process we only take the line flux from photometric data. 

There are 22 targets observed more than once. Four of them show no line but continuum 
and are excluded as LAEs. Two LAEs observed twice were confirmed with single frame's 
spectra, since their another spectra were affected by CCD gap or bright star contamination. 
The redshift values estimated from the line peak of the 16 multi-frame confirmed LAEs show 
consistence within 1 pixel . 

The four candidates show no line but continuum are all visible in the 2-d spectra of 
grism-300, while only two can be resolved in the grism-200 spectra. From table 2 we can 
see the first day observation with grism-200 show a relative low successful-fraction. This is 
mainly due to a lack of at least 2 ks exposure time compared with other frames. 

3.3. The stacked Spectrum 

We follow Wang et al. (2009) to stack the spectra of the confirmed Lya emitters in 
order to understand the nature of the line emitters. Excluding the spectra with very large 
background noise at the estimated line wavelengths, we stack 28 confirmed Lya emitters 
at z ~ 4.5 with grism 300 line/mm observation (see Figure [6]). We co-add the spectra in 
a variance-weighted method, and using a 2a clipping algorithm (one iteration) to remove 
artificial features(e.g., sky line residuals, ccd edges, etc) in the co-addition of each spectrum. 

The only visible line feature in the composite spectrum is the asymmetric Lya line 
(highlighted in the cutout of Figure |6]). We obtained wavelength-based asymmetry of ax = 
1.491q^3 fiux-based asymmetry of aj = 1.451;qq7, consistent with those expected from 
the asymmetric profile of high- redshift Lya emission line. The error bars on ax and af are 
estimated based on a large number of Monte-Carlo simulations, in which we added random 
noise to each data bin. 
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Fig. 4. — Plot of two measurements of the line asymmetry for 46 Lya emitters (see text 
for definitions of aA and a^.) The orange triangles and red stars are LAEs confirmed with 
10.2 ks exposure and 12.2 ks exposure with 200 line/mm grism observation, and the green 
diamonds and cyan squares are LAEs confirmed with 16.2 ks exposure and 14.4 ks exposure 
with 300 line/mm grism, respectively. 
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Fig. 5. — The Lya photometry hne flux compared with their spectroscopic Lya line flux. 
The orange and red stars are LAEs confirmed with 10.2 ks exposure and 12.2 ks exposure 
with 200 hne/mm grism observation, and the cyan and green diamonds are LAEs confirmed 
with 16.2 ks exposure and 14.4 ks exposure with 300 line/mm grism, respectively. 
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Fig. 6. — The co-add spectrum of selected 28 LAEs with grism 300 hne/mm spectra. The 
Lya hne is highhghted in the cutout, and the co-add spectrum at estimated position of CIV, 
Hell and Lyman series are also plotted. The errors of the co-add spectrum is plotted in red 
dotted line. 
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There is no possible He II A 1640 A or CIV A1549 A emission line near the expected 
wavelength (assuming the line width same to the Lja line width). At the expected wave- 
lengths of NV A 1240 A , the significance is about ~1 a. So neither line was detected, with 
upper limits <2a in all the cases. We determined the 2cr upper limits under the assumption 
that the undetected line has a physical width equal to the observed Lya width. The resulting 
limits on the line ratios are /(NVA1240)//(Ly-a) < 7.2%, /(CIVA1549)//(Ly-a) < 9.2% 
and /(HeIIA1640)//(Ly-a) < 6.5%. Our constrain on CIV is weaker than the previous work 
by Dawson et al. (<8%; Dawson et al. 2004) based on 11 LAEs spectra with Keck and 
Wang et al. (<3.7%; Wang et al. 2009) based on 110 LAEs spectra with IMACS, while the 
constrain on Hell is stronger than Dawson et al. (<13%) and comparable to Wang et al. 
(<7.4%). 



4. DISCUSSION 

4.1. AGN fraction and Unobscured Star Formation Rate 

The X-ray analysis on this sample has been done by Zheng et al. (2010) with the old 
2Ms CDFS data. We found one detection in ECDFS region, which has been spectroscopically 
confirmed by Treister et al. (2009) as a type 1 quasar at z = 4.48. X-ray stacking of the 
remaining sources shows a marginal detection (2.4cr), while spectroscopic data show that 
about half of the signal was from one LEG and one possible LAE (NB673-27). Excluding 
the two, we estimated a constraint (3cr) on the average X-ray luminosity of L2-skeV,rest < 
2.4x10^2 grg s"^ (Zheng et al. 2010). With the 4 Ms CDF-S data, we can decrease the 
average X-ray luminosity to L2-skeV,rest < 1.6x10^^ erg s~^. With the ratio of Lya to X-ray 
fiuxes for typical AGNs, we can estimate that fewer than 2.1% (4.2%) of our LAEs could 
be high redshift type 1 (type 2) AGNs, and those hidden AGNs likely show low rest frame 
equivalent widths. 

We also check the AGN properties (e.g, the broad emission lines or high ionized emission 
line such as CIV) in the spectra of our confirmed LAEs. No AGN lines have been identified 
in this work (Note that we did not target the 5 brightest LAEs including the quasar). We 
follow Dawson et al. (2007) and Wang et al. (2009) to use the line ratio of CIV and Lya from 
the composite optical spectra to constrain the upper limit of AGN contamination fraction. 
Wang et al. (2009) showed that the CIV line fiux is invisible even with stacking 110 Lya 
emitters at z ~ 4.5. They can only give a upper limit of the CIV to Lya line ratio of < 
3.7% (2cr). Compared with the typical ratio of CIV to Lya line in type II AGNs of 22% 
(Ferland & Osterbrock 1986), the AGN contribution to LAEs is <3.7/22 = 17%. In this 
work, we have presented the composite spectra of 28 Lya emitters at z~4.5 in Section 3.4. 
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The upper limit of our CIV to Lya line ratio is < 12.7% {2a), much weaker compared with 
the constraint of Wang et al. (2009). 

Deep radio data is also an independent method to check the existence of AGNs. In 
ECDFS field, there exists deep radio data taken with the VLA array and covered the entire 
ECDFS field to a typical sensitivity of 7.2 /iJy per 2.8" x 1.6" beam (Miller et al. 2008, 2011). 
The radio catalogue from Miller et al. was searched with a match radius of 3" with our LAE 
candidates and two objects were found. One object (J033127.2-274247, with a separation of 
<1") is the X-ray detected and spectroscopically confirmed z=4.48 quasar. The other one 
(NB673_21) with separations of 2.4" to radio detections, is contaminated by other nearby 
sources within 3" in their optical B-band image. So the X-ray detected z = 4.48 quasar also 
shows significant radio emission. The radio stacking of the remaining 111 LAE candidates 
reveal no-detection, and we get down to 0.67 uJy rms (~ 7.5/sqrt(110)). Since AGNs can 
be divided into radio-load and radio-quiet AGNs, and we don't know the relative number of 
radio-loud AGN and their characteristic radio luminosity, it is impossible to constrain the 
AGN fraction from the radio stacking results. 

Compared to AGNs, the star-formation activity in star-forming galaxies also contribute 
X-ray and radio emissions, and the X-ray and radio emission in star-forming galaxies are 
unobscured by dust (Zheng et al. 2011). Although the X-ray and radio radiation from 
galaxies are much fainter than that from AGNs, the number density of galaxies is significantly 
larger than that of AGNs. Xue et al. (2011) pointed out that in the 4Ms CDFS region, the 
X-ray radiation from galaxies begin to dominate the faint flux end in the log{N) — log{S) 
plot. If the X-ray radiation is all due to star formation (Ranalli et al. 2003), our X-ray 
average flux would correspond to a ?>a upper limit of unobscured star-formation rate (SFR) 
<320 M0 yr-^ If the radio fluxes are converted into star formation rate using the conversion 
rate of Yun et al. (2001), the rms sensitivity of the radio map corresponds to upper limit 
of SFR ~1700 M0/yr at z~4.5. The radio stacking of the remaining 111 LAE candidates 
can be translates to a 3-0" upper limit of SFR ~15OM0/yr at z~4.5. The average SFR from 
the Lya emission line (with the relation from Kennicutt 1998) is about 5 Mq yr~^. If we 
assume that the SFR from X-ray or radio is consistent with the intrinsic SFR, the ratio of 
SFRs from observed Lya and from upper limits of X-ray or radio can be used to constrain 
the Lya escape fraction, which is fesc{Lya) > 1.6-3.3% for z ~ 4.5 LAEs at 99% confidence 
level. 
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4.2. Population III stars 

Population 111 stars (or Pop 111 stars) are thought to have very strong Hell emission 
(Hell A 1640 line, the Balmer a transition of singly ionized helium). In this work, we get the 
upper limit of Hell-to-Lya ratio of < 6.5%, comparable with <7.4% of Wang et al. 2009. 
Our upper limit can be converted to ^o^^Que+IQu) < 1-46, where QHe+ and Qh are ionizing 
fluxes for He+ and H, respectively (Schaerer 2003). Wang et al. had pointed out that it is 
not easy to rule out the existence of Pop 111 stars in LAEs at z ~ 4.5, since the constraints 
only reach the maximum predicted values for metal-free populations with IMFs including 
very massive stars (Schaerer 2003), and the fraction of Pop 111 galaxies among LAE samples 
at z ~ 4.5 could be very small (e.g., Scannapieco et al. 2003). Recently, McLinden et al. 
(2010) reported the discovery of outflows in two LAEs at z~3.1, in which the Lya emission 
had a slightly higher redshift (outflow velocity ~ 300 km/s) than the rest-frame optical [0111] 
emission. In the co-add process above, we fix the Lya peak as the systemic redshift for each 
LAEs. If there exists outflows, there should be about <1.5 A blue-shift for the Hell line in 
the rest frame. So the existence of outflows will increase the uncertainty of Hell emission in 
our stacked spectrum. If we assume the line width as 5A in the rest frame, there is still no 
detection, and the upper limit of Hell-to-Lya ratio increases to < 8.4%. 



4.3. Ly« Luminosity Function 

We now present the Lya luminosity function (LP) based on our spectroscopically con- 
firmed LAE samples at z ~ 4.5 in ECDFS. We choose a modified version of the V/Vmax 
method in Dawson et al. (2007). Dawson et al. (2007) had considered two types of incom- 
pleteness in Lya LP measurement: the target incompleteness (i.e., not all candidates were 
targeted) and the spectroscopic sensitivity depth. The target incompleteness is the ratio of 
all candidates' number to targets' number at relative Lja flux bins (see figure [7]). We ig- 
nore the spectroscopic sensitivity incompleteness, as our spectroscopic depth are much more 
deeper than the photometry depth. We also take into account the narrowband incomplete- 
ness in our Lya LP calculation. The narrowband incompleteness was estimated following 
Hibon et al. (2010): we added 200 artificial star-like objects (LAEs are point sources in these 
NB images) per bin of 0.1 mag on the NB673 image, then ran SExtractor on this image for 
object detection, and got a direct measure of the narrowband incompleteness by counting 
the number of artificial stars detected in each magnitude bin. 

The comoving volume Vmax for each confirmed LAE was calculated, here Vmax was the 
volume where the source could be selected by our survey. The survey area is 34' x 33' from z 
= 4.44 to z = 4.56 (NB665 and NB673 filters) and 160 square arcmins from z = 4.36 to z = 
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-16.6 -16.4 -16.2 

Logio(fLya ergs s"' cm"^) 



-16.0 



Fig. 7. — Distribution of the Lya fluxes for all candidate LAEs (empty histogram), and 
those targeted (light grey histogram) and spectroscopically confirmed (dark grey histogram) 
LAEs. The dash-dotted line shows the spectroscopic successful fraction (i.e., the ratio of the 
confirmed number to the targeted number per flux bin). 
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Fig. 8. — Differential Lya Luminosity Function for our z = 4.5 sample, computed using the 
V jVmax method. The error bars are the 1 a statistical uncertainties given by the root variance 
of the 1/Knax- The black solid histogram gives the number of individual sources contributing 
to each luminosity bin. The dotted curve shows the best-fitting model for Wang et al. 2009 
with a fixed faint-end slope a = -1.5. The yellow, green and blue points are presented for 
narrowband NB656, NB665 and NB673 data, respectively. The vertical color lines show the 
detection limits at corresponding narrowband. 
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□ Z~4.5 W09 O Z~5.7 O08 
■ Z~4.5 DOT # Z~5.7 S06 




42.6 42.8 43.0 43.2 43.4 

logL*y„ [ergs/s] 

Fig. 9. — Contour plot of the fitting parameters L* and $* on tfie differential Lya Luminosity 
Function for our sample. Here red, blue and green contours are fitting results on our all, 
NB673 only, and NB665 only samples, with confidence level of 68% and 90% (A^^ = 2.3 and 
4.6). We compare L* and $* from various surveys at different redshifts. W09: Wang et al. 
(2009); D07: Dawson et al. (2007); O08: Ouchi et al. (2008); S06: Shimasaku et al. (2006); 
Kll: Kashikawa et al. (2011); Hll: Hu et al. 2011; OlO: Ouchi et al. 2010. Best-fitting 
parameters obtaned by fixing a = -1.5 are taken for all surveys except for D07 (a = -1.6). 
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4.43. The target incompleteness and narrowband incompleteness were taken into account. 

The derived Lya luminosity function of LAEs at z~ 4.5 in ECDFS is plotted in Figure 
El Following Malhotra & Rhoads (2004), we fit the luminosity functin with a Schechter 
function 



In fitting, we choose to fix the faint-end slope a = -1.5 since our data are not deep enough 
to constrain the Lya LF at the faint end. The best-fit parameters are log(L*) = 42.86 ±0.10 
and log($*) = -3.40±0.14 {x^/dof = 13.2/7). When applying the same method of the Lya 
LF measurements on the two narrowband images, NB665 and NB673, respectively, the best- 
fit parameters are log(L*) = 42.80 ±0.21 and log($*) = -3.42 ±0.33 {x^/dof = 0.6/3) for 
NB665 data, and log(L*) = 42.87 ±0.11 and log($*) = -3.21 ±0.15 {x^/dof = 30.3/7) for 
NB673 data. We plot the contour of the fitted L* and $* in Figure [9l which show some 
kind of coherence. The contour plot of fitting parameters on NB665 and NB673 differ at a 
confidence level P > 90%, implying a > 2a variance on L* and $*. Compared to previous 
work, the best-fitted L* and $* are consistent with the range of values found in various 
surveys at 3 < z < 6.5 (see Figure [9]) 

Ouchi et al. (2010) had found a ~2 a decrease on L* and $* from Lya LF at z = 5.7 
to z = 6.6 in the 1 deg^ Subaru Deep Field (SDF). They thought that their field is large 
enough to ignore the cosmic variance, so the ~2 a difference is probable due to the neutral 
hydrogen during the universe reionization. We found that their L* and $* from Lya LF at 
z= 5.7 and z = 6.5 are located in the 1-a contour area of our L* and $* from NB673 and 
NB665, respectively. And their comoving volume of 8x10^ Mpc~^ and 1x10^ Mpc~^ at z 
= 5.7 and z = 6.5 in the SDF, respectively, are about 3~4 times larger than our volume of 
2.3x10^ Mpc~^ per narrowband image at z = 4.5. So the 2 cr difference of Lya LF at z = 
5.7 and z = 6.5 may be also due to the cosmic variance. 



We present a spectroscopically confirmed LAEs at z ~ 4.5 in the ECDFS region. This 
sample is from two contiguous narrowband images (NB665 and NB673), and a much smaller 
region of a shallower narrowband image (NB656). The main scientific results are summarized 
as below: 

• We identify 46 z ~ 4.5 LAEs, 5 LBGs and 1 [OIII] emitters from our targeted 68 
LAE candidates. The Lya conrmation fraction is about 70%-80%, and the contamina- 
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5. 



CONCLUSION 



- 21 - 



tion fraction is about 14%-19%. All targets with fiya > 3.7x10 ergs cm ^s ^ are 
confirmed. 

• We do not find any CIV or Hell lines over 2-cr significance in the spectra of the 
confirmed LAEs, even with stacking methods. 3 

• There is one quasar previously confirmed as type 1 AGN at z = 4.48 in our sample, 
which is also detected with the deep X-ray and radio observations. The remaining 
LAEs are not detected in X-ray or radio even with stacking methods, and the stacked 
upper limits in radio and X-ray can be converted to average SFR as SFRija^io < 150 
M0 and SFRx < 200 M©, respectively. 

• The upper limit of Hell to Lya line ratio in our coadded spectrum can be converted 
to Qhc+IQe < -1.46. However, we can not exclude the existence of Pop HI stars, and 
the existence of outflows will increase the uncertainty in determining the real position 
of Hell line. 

• The Lya luminosity functions in our sample are consistent with the previous surveys at 
3 < ^ < 6.5, as lg(L*) = 42.86 ± 0.10 and lg(<I>*) = -3.40 ± 0.14 with x^ldof = 13.2/7. 
However, the Lya LP differs at > 90% confidence level in the individual narrowband 
images, implying the existence of cosmic variance. This difference is comparable to the 
Lya LP evolution at z = 5.7 and 6.5. 



We would like to thank the support of NSP grant AST-0808165 and NOAO TSIP 
program. The work of JXW is supported by National Basic Research Program of China 
(973 program. Grant No. 2007CB815404), and Chinese National Science Foundation (Grant 
No. 10825312, 10773010). 



Table 3. Spectroscopically Confirmed Lya Emitters at z~ 4.5 in the CDF-S field 
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^ The Lya line flux is calculated from the narrowband and broadband photometry in the unit of 10 
erg cm~2. 

LAEs with non-detection in the R-band are marked with EW > 10000 A. 

LAEs confirmed from their spectra are assigned grades, here 1 means very good quality, 2 means not 



very significant, and 3 means low s/n. 

The Lya line widths are directly measured from their spectra and not corrected for the instrument 
profiles. 
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